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ABSTRACT 

If cooled-down neutron stars have a thin atomic crystahine-iron crust, they must diffract X- 
rays of appropriate wavelength. If the diffracted beam is to be visible from Earth, the illuminating 
source must be very intense and near the reflecting star. An example is a binary system composed 
of two neutron stars in close orbit, one of them inert, the other an X-ray pulsai0. The observable 
to be searched for is a secondary peak added (quasi-) periodically to the main X-ray pulse. The 
distinguishing feature of this secondary peak is that it appears at wavelengths related by simple 
integer numbers. A, A/2, A/3... X/n because of Bragg's diffraction law. 

Subject headings: Bragg difraction, neutron stars, crystalline iron, Bragg peaks 



Introduction 



Cu rrent theories of Neutron Stars (IBombaci 
20071) imply that as pressure builds up towards 
the interior of the star there are successive phase 



transitions from an iron crust ( Chamel fc Haensel 



2008t Ide Young 111991 ). to a nuclear medium with 



high neutron density, to a neutron Fermi liquid, to 
more exotic forms of matter. It has proven difficult 
to make empirical progress on the star's compo- 
sition. For ex ample, very informa t ion-rich equa 



tions of state (Huber et al 



1994 iKlahn I [iooeh 

have to be tested with few numbers (mass, size, 
pulsar period and time dependence, etc.) 

That the outer-most layer of the star contains 
iron is know n from the cha r acteri stic absorption 
lines of Fe ( Cottam et al. ] I2OO2I) . In recently 
formed stars or in stars that are heated by ac- 
cretion, a hot iron "ocean" is likely to cover the 
surface ( Haensel et al. I2OO6I) . However, if the star 
is not accreting material, it cools rapidly; the star 
stops emitting X-rays afte r about a million years, 
(jChamel fc Haensel |[2008l) . At about 2 x lO'^ years 
(a small time in galactic scale) the star is be- 
lieve d to have reached 1000 K simply by radia- 
tion (lYakovlev fc Pethick1l2004t llmshennik et~ar 



20021 ). a temperature which is well below the so- 
lidification of iron. 

The hypothesis that this iron layer is in crys- 
talline form has only been indirectly tested by ob- 



servations of the initial cooling in quasi-persistent 
soft X-ray transients. These seem to be consistent 
with the neutron star crust having the structure of 
a perfect crystal, while models ba sed on an amor 



phou s crust cannot fit the data (jShternin et al 



20071) . 



We here propose a direct test that might be per- 
formed as the catalog of X-ray sources expands, 
requiring the existence of a binary system com- 
posed of one such cooled-down neutron star with 
certain minimal assumptions that will be spelled 
out shortly, and an X-ray emitting companion. 

The canonical way of ascertaining the crys- 
talline structure in a laboratory material is by 
exposing it to X-rays and studying the resulting 
diffraction pattern (it is not possible to directly 
observe neutron diffraction from Earth due to the 
neutron's short lifetime of about 15min). This 
textbook method works well even for fragmented 
crystals, for example a dust sample made of mi- 
crocrystals can be made to produce a diffraction 
pattern. 

The concept is the same in an astronomical con- 
text: the "sample" is the neutron star whose crys- 
talline crust is being examined, and one needs a 
beam powerful enough to detect it, after reflec- 
tion on the star, in Earth or in a orbiting X-ray 
telescope. This stringent condition requires a very 
intense beam to be focused on the star and then 
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reflected onto our direction. 

A promising and simple system that offers a 
potential opportunity to explore this diffraction 
is a binary system composed of two neutron stars. 
We will assume one of them, nicknamed ^^Pharus" , 
to be an active X-ray pulsar directly visible from 
Earth. Many pulsars ar e known to emit in X-rays 



(jCampana et al. Its period can be taken 



of order 1-second, as typical for the larger part of 
the normal pulsar population. Its companion in 
the binary system, that we will name Reflector^^ 
is assumed to be another compact object, indeed 
a neutron star, but this one dark and inert, held 
by its neutron degeneracy pressure but emitting 
no significant radiation, nor possessing a strong 
magnetic field that may alter the refiection of the 
beam from Pharus. If the orbital conditions are 
right, then one may observe the direct pulses from 
Pharus, but also at least one diffracted beam from 
Refiector, depending on the frequency of X-ray ob- 
servation. 

A further complication that may arise is if the 
neutron star possesses an atmosphere of ionized 
gas above the thin iron crust that we would like 
to see detected. This po ssibility has been analyze d 



with quite some care in ( Suleimanov et al. 2009() 



In fact, part of the X-rays in a magnetized star 
are scattered or absorbed by this atmosphere. It 
is not clear whether this is a general phenomenon 
and applies also to a non-magnetic star, but it 
nonetheless needs to be kept in mind. A double 
(binary) pulsar would be more difficult to ana- 
lyze since in addition to the main and reflected 
pulses there would be a second main pulse set and 
maybe even another set of secondaries. Moreover 
the strong magnetic field of the mirror star could 
distort the reflecting crystal. 

Of course many other forms of intense X-ray 
emission are possible, such as bursts from accre- 
tion processes in X-ray binaries or from black - 
hole or supernovae ejecta (jBorkowski et al. II2OO7I) . 
These should also provide the equivalent of Bragg 
peaks (but without the characteristic main pulse 
of Pharus) . An advantage of a binary neutron-star 
system is that the main pulse provides a clock- 
time T'^P'" that is useful in the data analysis. 



2. Orbital conditions for observability of 
Reflector 

We study a simple orbital case, in which the two 
stars have equal mass and are therefore equidistant 
to the center of mass. This is an extreme situa- 
tion and, on a case by case basis, one could study 
more general configurations should data become 
available. Known pulsars forming part of dou- 
ble neutron stars have masses that cluster around 
1.35 sol ar masses with a spread of order 0.04 solar 



masses (jThorsett fc Chakrabartvl Il999l ). There 



fore this case is close to what one might expect to 
find when looking for a double neutron-star sys- 
tem. 

Pharus presents, as typical for pulsars, a strong 
magnetic field precessing around its spin axis. X- 
rays are emitted along the direction of the mag- 
netic field lines streaming out the magnetic pole; 
thus, the radiation beam, which is observed as a 
time series of hashes as it rotates, swipes a conic 
surface. The visual to Earth, PT in figure[Tl has to 
be contained in this surface, since, after all, Pharus 
is assumed to be observed as a conventional X-ray 
pulsar. 

Further, we will assume that the plane of the 
orbit of Reflector (taken circular) contains the vi- 
sual to Earth, so that the orbit is seen on edge. 
These are not necessary conditions, and a com- 
puter programme can be written to incorporate 
more general cases if putative data becomes avail- 
able. The situation is depicted in figure [1] 

The beam cone intersects the orbit of Refiector 
in at least one point, and its angular position along 
the orbit will depend on the opening of Pharus's 
cone (angular width), controlled by its magnetic 
anomaly or inclination angle to the rotation axis. 

If Refiector happens to be at the orbit-beam 
cone crossing when Pharus illuminates it, the re- 
flected radiation might reach Earth at an instant 
when the pulsar is, in theory, still. 

In our simple computer simulation, the inten- 
sity of the radiation of Pharus for a given X-ray 
wavelength will be taken as the unity at maxi- 
mum, with a Gaussian profile (this plays no role 
as long as it is a simple function). The intensity 
of the diffracted beam will naturally be smaller. 
This will be partly due to absorption in Reflector, 
but primarily because the illumination of its sur- 
face will decrease with the square of the distance 
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to Pharus, 




V 

T 



Fig. 1. — T denotes the terrestrial observer; P 
the pulsating star Pharus; R the reflecting star 
Reflector; if R in its orbital motion crosses the cone 
illuminated by P's X-ray beam at a time where 
the crossing is illuminated, then Bragg diffraction 
might occur, visible at a wavelength meeting the 
Bragg condition for the reflecting angle 9. 



i?2 



(1) 



where r is the effective emitting radius of Pharus 
(of order 100 km for a typical neutron-star pulsar 
of radius 10 km) and R the radius of the orbit of 
Reflector. 

To keep the signal of Reflector within one or 
few percent of the intensity of Pharus, Reflector 
has to orbit in close range to Pharus. However, to 
have a clearer line-shape for the signal it is conve- 
nient to have the orbital period of Reflector, T°'^^^^ 
to be of the same order as, (but larger than) the 
spinning period of Pharus, T'^p'" ~ Is, which de- 
mands, because of Kepler's third law. 



T 



' orb it 



y^{Ml+M2jd 



, if M2 = Ml 



larger orbital radii. Given the tension between 
the two requirements, we compromise to an or- 
bit of order 10000 km (or a period T°''''** of order 
15s), which is certainly small enough for relativis- 
tic corrections to the orbit to play a role, but we 
will ignore these. (Note that the double neutron 
star binary with the shortest orbital period is, a s 
of 2004, 2.45 hours in J0737-3030 (jLvne et al. 



20041 ) ■ so that this case study is low compared 
with current observations of binary systems, but 
not impossible). 

To check that relativistic corrections are not a 
major source of uncertainty, we take the relevant 
equations of General Relativity (Peters 1964) to 
estimate in turn the power radiated by gravita- 
tional waves, the rate of decrease of the semi ma- 
jor axis, the rate of decrease of the orbital period, 
and the total time for collapse of the binary, yield- 
ing, for typical values of the star's mass and the 
period/orbital radius just quoted. 



P=^^^Mi^f^^L±^^10-Watt (2) 



5 c5i?5 
dR _ 64 2G^M^ 



0.01 - O.lm/s (3) 

: (4) 



dt 

96 Ml M2 (Ml +M2) . 47r2 4 

5 g5 (T^orbit ) I ^G{Mi + M2) 

1 — 10/is/year 



3 



Thus, relativistic effects may be neglected if tlie 
observation time of pulsar profiles is very short 
compared with the time to collapse just estimated. 

The extraction of the reflection angle for a 
given binary system might employ the variation of 
Pharus' velocity along the line of sight in its orbital 
motion, causing a first-order kinematic Doppler 
shift, affecting the spectral shape of the main peak. 
Analyzing the red/blue shift of the peak over sev- 
eral periods one can determine from the data to 
what angular position along the orbit does a given 
main peak of Pharus correspond. Then one just 
needs to note that the position of Refiector a short 
time later takes a correction 



_ g|Doppler _|_ '2nAt 
^period 

3. Bragg diffraction at Reflector 



The principle of Bragg diffraction is illustrated 
in figure [2j A crystalline material acts as a difrac- 




sin0 



^collapse ^ 5 _ lOyears (5) with d the spacing between crystal planes, A the 

radiation wavelength, 9 the reflection angle taken 
from the vertical, and n an integer number. Be- 
cause the sine is smaller than one and n greater, 
this condition can only be satisfied for wavelengths 
smaller than twice the crystal spacing. 

It is known that atomic iron crystals present 
a body-centered cubic structure with spacing d = 
0.3nm called ferrite at temperatures above 1665 K 
or below 1184 K. Austenite, a crystallization 
in a face-centered cubic lattice is stable between 
these two temperatures (jKawasaki steel I |199Z) . 
These cubic crystallizations are already quite op- 
timal and therefore we can assume that the out- 
ermost layer of a neutron star, if crystalline iron, 
will employ them depending on the temperature. 
Therefore, if the spacing d is not altered in neutron 

(6) stars, the characteristic wavelengths that need to 
be scanned are A < O.Gnm. 

That the lattice spacing remains, at the neu- 
tron star's surface, similar to that under earthly 
conditions requires closer examination. First, we 
show from a macroscopic point of view that grav- 
itational forces at the surface of Reflector only 
moderately deform the crystal. Indeed, the in- 
wards gravitational acceleration at the surface is 
of order g = lO^^m/s^. This causes a force per 
unit area F/Aq — gph to depth h inside the 
crust (the density of iron being p — 7.87 g/cm^). 
Given the Young modulus of iron, Y — 196.5GPa, 
the distortion of the lattice expected is IS.L/L — 
F/{Y Aq), which grows linearly with depth. The 
penetration depth of X-rays depends on their fre- 
quency and is typically 2 — 15/im in crystalline iron 
(|Randle k Engler 1 12000| ) (this amounts to a num- 
ber of about 10^ — 10^ planes necessary to totally 
diffract the beam). AL/L remains signiflcantly 
smaller than 1 for all layers within this penetra- 
tion depth. 

Another argument is to consider that char- 
acteristic microscopic atomic forces are of order 
2.5 X 10~® N and these, based in Coulomb inter- 
actions between electrons and iron nuclei, are also 
much stronger than gravitational forces exerted on 
the atomic mass, of order 3 x 10~^^ N. Thus, for 
the purposes of X-ray diffraction, crystalline iron 
at the star's surface is similar to that on Earth, 
with somewhat broadened lines. 

However, the radiation wavelength that we ob- 

(7) serve is red-shifted respect to the radiation re- 



Fig. 2. — If the radiation's wavelength is com- 
parable to the crystal spacing, the crystal acts 
as a diffraction grating, reflecting different wave- 
lengths by different angles. The phenomenon is 
understood as constructive interference, with the 
rays (black arrows) having all equal phase at the 
wave front (transverse red lines). The additional 
path travelled by rays penetrating deeper into the 
crystal needs to be an integer number of wave- 
lengths for the phase of these rays to match those 
reflected at the surface. Bragg's law follows. 

tion grating. A beam of a given wavelength is 
reflected, by constructive interference, only when 
Bragg's law is met 



nX 
2d 
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fleeted at the star's surface, due to the gravita- 
tional potential in General Relativity. This eff'ect 
brings about a correction necessary when extract- 
ing the crystal spacing from Bragg's law, since the 
A appearing in that formula is taken at the star's 
surface. One finds 



A 



Earth 



A(l + z) 



A 



v/1 - 2MG/RC 



(8) 



in terms of the star radius and mass. (This cor- 
rection factor is about 1.3 for a typical neutron 
star). 

Equation ([7]) would then be corrected due to 
the gravitational red-shift to 



2d{l + z) 



(9) 



This formula can also be corrected for kinematic 
Doppler shift due to the orbital motion of Reflector 
as neeed be, simply substituting the factor {1 + z). 

If a diffraction pattern was suspected, one could 
conceivably identify the form of crystallization 
(whether ferrite or austenite) by the lines present. 
This is possible because the fee and bee crystal 
have different symmetry planes. In terms of the 
Miller indices /i. A:, / and the unit-cell size a, the 
spacing between those planes is given by 

(10) 



/l2 + B + P 

and combining with Eq.® 



sin0 



2a(l + z) 



(11) 



For a bee lattice only those combinations where 
h + k + 1 is an even integer occur in the diffraction 
pattern. An fee lattice diffracts instead for h, k, I 
all odd or all even. Table [T] shows the first few 
diffraction lines expected. Given are the Miller in- 
dices, the distance between diffracting planes with 
given Miller indices in lattice units, and the mul- 
tiplicity of planes yielding the same distance (and 
thus diffracting at the same wavelength), the in- 
tensity of the diffracted beam being proportional 
to this multiplicity. Additional factors influenc- 
ing the intensity of a given line are the Lorentz- 
polarization factor and the atomic scatteri ng fac- 
tor for Fe, given in (|Cullitv k Stock Il200l[ ). 



Table 1: Inverse of the interplane spacing relative 
to the unit-cell size a/d and multiplicity factor P 
(how many combinations of Miller indices lead to 
the same d) for the fee (left) and bec(right) lat- 
tices. 
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Fig. 3. — The spectrum of the binary system for the case of equal masses. Here the pulse for several 
revolution periods of the system (about 15 seconds each). The repetition rate of the reflected pulse depends 
on the orbital period of Reflector and the spinning (pulsar) period of Pharus, since Reflector has to cross 
Pharus's beam cone at a moment when this is illuminated. The intensity of Reflector relative to Pharus 
(taken as 1 at its maximum) depends on the albedo of Reflector as well as on its distance to Pharus. It 
could be much larger if the beam of Pharus was very well collimated. Note that not all pulses have the same 
height due to the different illumination of Reflector in different crossings of the radiation cone of Pharus. 
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4. Line shape of the pulsar 

An example of an observable line-shape for the 
binary X-ray pulsar is plotted in figure [31 

One sees that, in addition to the pulse of 
Pharus, a smaller signal due to Reflector appears. 
The intensity of this pulse depends on the distance 
between Pharus and Reflector, but Kepler's third 
law links this to the orbital period, and thus the 
relative intensity of the secondary pulse is not in- 
dependent of its repetition rate. 

This periodic reapparition of Reflector depends 
on the orbital coincidence that it crosses the cone 
subtended by Pharus's beam when it is irradiated. 
This happens when, for two integers h and I2, 

At = hT°'^'^ = hT'P''' . (12) 

This rational number relation does not need to be 
exact, in view of the time- width of Pharus's beam. 
We take some ad-hoc values for illustrational pur- 
poses that are Tortu — 14.89 sec and Tspin = 1 sec. 
The resulting line shape, where the secondary is 
clearly observed, is represented in figure |3l Note 
the secondary misses several periodic apparitions 
since the spin and orbital period are not perfectly 
matched (which we expect to be borne in reality). 
In figure [4] we represent the same data for shorter 
periods so one can discernt the Gaussian shape of 
the main pulse. 

What distinguishes Bragg diffraction from 
other astrophysical phenomena that produce com- 
plicated secondaries (i ndeed pulsars have very dif- 
ferent line-shapes t Graham- Smith 1 120031 ) is that 
the peak appears at very specific wavelengths. 
Since n in eq. (O is an integer, the secondary 
appears at the wavelength sequence 

AAA 
' 2' 3' 4'" 
and at other specific frequencies as dictated by the 
various planes in table [TJ 

In addition, for the degenerate case (orbit seen 
edge-on) that we have chosen, there is one addi- 
tional crossing of Pharus's beam cone with Re- 
flector's orbit. The reflection angles appearing in 
Bragg's law for these crossings satisfy, given that 
the binary system subtends a negligible angle as 
seen from Earth, 29i + 202 — tt. Then one can 
easily see that a second series of wavelengths 

y V y 

2 ' 3 ' 4 



with A ^ = (2(i)^ — A^ will also show a secondary 
peak (same for other sequences caused by addi- 
tional reflecting planes). 

Another very interesting degenerate case occurs 
when the magnetic axis of Pharus is perpendicular 
to its precession axis: then the X-ray beam swipes 
the entire orbit of Reflector, and yields the pos- 
sibility of detecting, at short X-ray wavelengths, 
several diffraction maxima. 
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Fig. 4. — Same as in figure [3] but detailing the 
combined pulse for short times. 



5. Conclusions and outlook 

In conclusion, we have attempted to develop 
an astronomical observable able to test the hy- 
pothesis of a crystalline iron crust in cold neutron 
stars. The search we suggest is for a secondary 
pulse in an X-ray pulsar indicating a binary sys- 
tem. The tell-tale of crystalline structure is the ap- 
pearance of this secondary for a specific sequence 
of wavelengths A, ^, j... Several such fre- 
quency sequences may appear due to the various 
crystal planes and to the possibility of observing 
Reflector at more than one angle along the orbit. 

The phenomenon is quite the opposite of an 
eclipsing binary. Instead, here the companion star 
is dark most of the time and only brightens when 
Bragg's diffraction condition is satisfied. 

To extract the crystal spacing d from data one 
could follow the (simplified) steps 

• Detect a potential binary system with the 
proposed pattern 
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• Measure the period of Pharus's main pulse, 

^spin 

• Perform a simple rational number analysis to 
obtain an approximate for Reflector. 

• Obtain the angle by a Dopplcr measure- 
ment or possibly a direct orbital reconstruc- 
tion. 

• Given Aq, the largest wavelength where Re- 
flector is seen, after any corrections needed, 
use Bragg's law to obtain d. 

The rate of coalescence of two-neutron star bi- 
nary systems has been estimated in several works 
iKim 2003; Regimbau 2005.). given its interest 
for the gravitational wave detection programme at 
LIGO and VIRGO, as well as possible upgrades of 
these machines. 

Estimates for double neutron star mergers 
inside our galaxy are (in order of magnitude) 
10~^ yr'"^ . To obtain a good rate that allows 
a potential detection in, for example, a decade of 
operation, one needs to be able to detect the X-ray 
source to a distance of 30-40 Mpc (more than a 
hundred times the galactic size) to include many 
other star formations of the local group. This is 
the goal for the gravitational wave detectors. 

The detection of X-ray diffraction in a binary 
system would point out to the tightness of the or- 
bit and the temporal approaching of its final colli- 
sion. Therefore it would be interesting to explore, 
in future work, the possibility of employing X-ray 
diffraction as an early-warning system for those 
detectors. 

X-ray diffraction might also occur at the surface 
of a cold white-dwarf acting as Reflector. A differ- 
ence between such stars and neutron stars is that 
nuclear fusion chains are generally thought not to 
be complete in the progenitor star, one does not 
expect to find iron but a mixture of lighter ele- 
ments such as oxygen, neon, silicon, etc. It would 
surely be interesting and complicated to examine 
any data that may become available. Immediately 
one notices that the observation of the iron crust 
on a cold neutron star is an additional test of the 
efficiency of fusion in its very massive progenitor 
(believed to be completely efficient). 

The authors thank Krysty Dyer and Katja 
Waidelich for a careful reading of the manuscript. 
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